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r ^ , Abstract. We examine the properties of central dust in nearby quiescent and active early-type galaxies. The 

active galaxies are low-power radio galaxies with Fanaroff & Riley Type I or I/II radio jets. We focus on (a) the 
£C) ' comparison of the dust distributions in the active and quiescent galaxy samples and (b) the relation between 

the radio jet and dust orientations. Our main observational conclusions are: (i) in line with previous studies, 
the dust detection rate is higher in radio-jet galaxies than in non radio-jet galaxies; (ii) radio galaxies contain a 
higher fraction of regular dust 'ellipses' compared to quiescent galaxies which contain more often irregular dust 
distributions, (iii) the morphology, size and orientation of dust ellipses and lanes in quiescent early-types and 
active early-types with kpc-scale radio jets is very similar; (iv) dust ellipses are aligned with the major axis of the 
galaxy, dust lanes do not show a preferred alignment except for large (>kpc) dust lanes which are aligned with 
the minor axis of the galaxy and (v) as projected on the sky, jets do not show a preferred orientation relative to 
' the galaxy major axis (and hence dust ellipses), but jets are preferentially perpendicular to dust lanes. 

We show that the dust ellipses are consistent with being nearly circular thin disks viewed at random viewing 
angles. The lanes are likely warped dust structures, which may be in the process of settling down to become 
regular disks or are being perturbed by a non-gravitational force. We use the observed dust-jet orientations to 
constrain the three-dimensional angle #dj between jet and dust. For dust-lane galaxies, the jet is approximately 
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perpendicular to the dust structure, while for dust-ellipse galaxies there is a much wider distribution of #d.i- 
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then the jet orientation apparently is roughly aligned with the angular momentum of the dust before it settles. If 



We discuss two scenarios that could explain the dust/jet/galaxy orientation dichotomy. If lanes are indeed settling, 



lanes are perturbed by a jet-related force, it appears that it causes the dust to move out of its equilibrium plane 
in the galaxy into a plane which is perpendicular to the jet. 
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■ Key words. Galaxies: active - Galaxies: elliptical and to power an active galactic nucleus (AGN). In the nearby 
lenticular, cD - Galaxies: nuclei - Galaxies: jets - ISM: Universe, many, perhaps the majority, of nearby early- 
dust, extinction. type galaxies display some sort of nuclear activity, most 

often as a modest AGN (e.g., Ho et al. 1997; Kauffmann 
et al. 12008(1 . Comparing the demography of the fuel reser- 
voirs in nearby active and quiescent galaxies can shed light 



1. Introduction 



Many early-type galaxies harbour dust and warm (i.e., on the triggering and/or feeding mechanism of AGN. In 

T~ 10 4 ~ 5 K) gas in their central regions (e.g., Sadler a small fraction of the active galaxies the AGN is accom- 

& Gerhard 1X3551 Goudfrooij et al. IT331 van Dokkum Pame d by kpc-scale radio jets (e.g., Condon & Broderick 

& Franx [RKty . Evidence is mounting that all early- EH- Here we focus on those radio galaxies. The orienta- 

type galaxies harbour a central supermassive black hole tion and kinematics of the dust and gas form a tracer of the 

(BH) (e.g., Kormendy & Gebhardt IMTfl and references stellar potential in addition to the st ellar p hotometry and 

therein). The dust and gas form potential fuel for the BH kinematics (e.g., Merritt & de Zeeuw 11983J). This requires 

that the interstellar material has settled in the gravita- 

Send offprint requests to: G.A. Verdoes Kleijn tional potential and is unperturbed by collisional forces. 

* Based on observations with the NASA/ESA Hubble Space The observed distri bution of relative orientations of the 



gas and dust, the jet, and the stellar potential, can con- 
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of dust and jet and perhaps the jet formation mechanism 
in radio galaxies. 

Many ground-based studies have surveyed the global, 
kpc-scale, dust content of nearby (distance < 100 Mpc) 
quiescent and/or active early-type galaxies (e.g., Sadler & 
Gerhard 113551 Veron-Cetty & Veron Knapp et al. 

H3B3I Roberts et al. ITMT1 Goudfrooij et al. 115911) . With 
the Hubble Space Telescope (HST) it has become pos- 
sible to detect routinely central dust distributions, also 
those with sizes below 1 kpc, in nearby galaxies. Using 
HST/WFPC2 optical broad-band photometry, dust is de- 
tected in ~ 50% of the nearby early-type galaxies (van 
Dokkum & Franx [T555I Tomita et al. I50UU1 Tran et al. 
200 lj. This detection rate increases to ~ 90% for nearby 
early-type galaxies with radio jets (van Dokkum & Franx 
119951 Verdoes Kleijn et al. H999f) . The increased detection 
rate in radio galaxies suggests a causal connection between 
extended central ISM distributions and the on-set of nu- 
clear activity and radio jet formation. 

Early studies, using ground-based imaging, showed 
that jets in radio galaxies are roughly perpendicular in the 
plane of the sky to kpc-scale dust structures (e.g., Kotanyi 
k EkersEZi Mollenhoff, Hummel & Bender H^jl . No re- 
lation was found between galaxy and jet orientation (e.g., 
Battistini et aL lTMH Birkinshaw & Davies|l985 bansom 
et al. H987[l apart from a tendency to avoid large *J 70° an- 
gles between galaxy minor axis and radio jet. In studies of 
the sub kpc-scale dust, based on HST imaging, the perpen- 
dicularity of jet and dust in the plane of the sky was also 
found for radio galaxies (e.g., van Dokkum & Franx 1995 
de Koff et al.EHOOl de Ruiter et al.|MEl>- Furthermore, 
evidence was reported that also the intrinsic, i.e., three- 
dimensional, orientation of radio jets is roughly perpen- 
dicular to the dust, using samples of radio galaxies which 
contain dust disks (Capetti & Celotti [1.999 , Sparks et al. 
2000). In contrast, Schmitt et al. lj2"0T)2"|l found m a sam- 
ple of 20 radio galaxies with regular dust disks that the 
jets are not roughly perpendicular to the disks in three- 
dimensional space. 

We have studied a sample of nearby Fanaroff & Riley 
ljT9T2jl Type I radio galaxies, the 'UGC FR-I sample' (Xu 
et al. ESnnni Verdoes Kleijn et al. IT5591 IM)2l Noel-Storr 
et al. |2T)U5J> . Verdoes Kleijn et al. ljll?9"9"jl describe a di- 
chotomy in apparent dust morphology of the UGC FR-I 
sample galaxies, which either contain regular dust struc- 
tures with an elliptical appearance ('ellipses'), or filamen- 
tary structures ('lanes'). Ellipses turn out to be aligned 
with the galaxy major axis while the lanes show no rela- 
tion to galaxy orientation. Furthermore, most lanes (but 
only some ellipses) are perpendicular to jets in the sky- 
plane. Clearly, the apparently conflicting results on the 
dust-jet orientation mentioned above might be connected 
to this relation between dust orientation and morphol- 
ogy. We explore this relation further by inferring the in- 
trinsic three-dimensional relative orientations of jet, dust 
and stellar content of the host galaxy. With the results of 
this analysis we hope to constrain the causal connections 
which lead to the observed orientations between the fol- 



lowing three actors: (i) the gravitational potential, (ii) the 
central ISM and (hi) the radio jet. We focus on the fol- 
lowing questions. Are the properties of the fuel reservoirs 
different in radio galaxies and non-radio galaxies? Is the 
dust orientation governed by the gravitational potential? 
Do jet/AGN related forces affect the dust orientation? 

This paper is organized as follows. In Section|21we de- 
scribe the selection and properties of the radio and non- 
radio galaxy samples. In Sections 01 and 01 we compare the 
dust properties in radio and non-radio galaxies as derived 
from HST/WFPC2 imaging. In SectionElwe describe the 
relative orientations of jets, dust and the galaxy in the 
plane of the sky. In Sections and we determine the 
constraints placed by the projected orientations on the in- 
trinsic relative orientations in three dimensions for galax- 
ies with dust classified as ellipses. We consider the dust 
lane galaxies in Section |HJ In Section we discuss the pos- 
sible origins for the dichotomy in intrinsic orientation of 
dust and jets found in Sections and |S] We summarize 
our results in Section ITUl 

Throughout the paper we use a Hubble constant Hq= 
75 kms^Mpc- 1 . 

2. Samples and Data 

We use a complete sample of 21 radio (i.e., 'active') galax- 
ies, the UGC FR-I sample, and a sample of 52 non-radio 
(i.e., 'quiescent') galaxies, the UGC non-FR-I sample, to 
statistically compare dust distributions. We also consider 
a larger 'FR sample', which contains 47 FR-I and FR- 
I/II radio galaxies in total, to improve the constraints 
on the three-dimensional relative orientation of jets, dust 
and stellar hosts. We rely mainly on HST/WFPC2 imag- 
ing. The exceptions are 3C 402N, for which we use 
HST/FOC imaging, and NGC 5128 for which we employ 
VLT/VIMOS imaging which provides a full view of the 
kpc-scale dust distribution in this nearby galaxy. 

2.1. The Radio Galaxy Samples 

The UGC FR-I sample is a complete sample of nearby 
radio galaxies from the UGC catalogue (Nilson et al. 
HWSjl . It contains all 21 nearby (v < 7000 km s _1 ), ellipti- 
cal or SO galaxies in the declination range — 5° < 6 < 70° 
in the UGC catalogue (magnitude tub < 14'. n 6 and an- 
gular size p > lfO) that are extended radio-loud sources, 
defined as larger than 10" at 3<r on VLA A- Array maps 
and brighter than 150 mJy from single-dish flux-density 
measurements at 1400 MHz. These galaxies have jets 
classified as Fanaroff & Riley ifO??^ Type I (FR-I). HST 
imaging for the UGC FR-I sample is presented in Verdoes 
Kleijn et al. ifTMl I2UU2|| . Xu et al. <|2UUU|) and Noel-Storr 
et al. I|2003| report VLBA radio imaging and HST optical 
spectroscopy, respectively, for the sample. 

We searched the literature for additional galaxies with 
FR-I or FR-I/II radio sources which show dust on HST 
imaging data. This resulted in 26 galaxies from the 3CR 
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and B2 radio galaxy catalogues and one 4C galaxy. We re- 
fer to this extended sample as the FR sample, and list the 
general properties of the 47 FR-I galaxies in Tabled The 
HST imaging data for 3CR and B2 galaxies are described 
in de Koff et al. (|PM)jl . Martel et al. ijTMSjl . Capetti et 
al. {2000) and Sparks et al. l|2000) while the imaging for 
4C-03.43 is described in Schmitt et al. (|2TiU2"|) . 

2.2. The Non-Radio Galaxy Sample 

Figure^Jshows the location of the UGC FR-I sample mem- 
bers among the galaxies of the UGC in the plane of cen- 
tral stellar velocity dispersion versus host absolute mag- 
nitude. The data are taken from the LEDA catalogue 1 
except for the velocity dispersion of 3C 66B and 3C 449 
(B alcells et al. UWty and NGC 4335 (Verdoes Kleijn et 
al. I2UU2)1 . The location of the UGC FR-I galaxies con- 
firms the well-known fact that radio galaxies are bright 
early-type galaxies. The figure also shows that the FR-I 
galaxies are distributed fairly uniformly, relative to the 
general galaxy distribution, in the upper-right corner of 
large host magnitude and velocity dispersion. The mor- 
phologies of the UGC FR-I galaxies are not a random 
selection from the morphology distribution in this upper- 
right corner. Table [21 shows that the UGC FR-I galaxies 
are predominantly ellipticals and E/SO transition objects 
and rarely pure SOs. 

To carry out a statistical comparison of the dust dis- 
tributions in early-type galaxies with and without radio 
jets, we need samples of both classes of galaxies with 
similar global properties of their stellar host and similar 
HST imaging of the central regions so that we can per- 
form an identical dust analysis. Ideally, one would like 
to perform the dust analysis for a well-defined sample 
of quiescent galaxies in the UGC catalogue from which 
also the radio galaxies were selected. Tran et al. (2001) 
analysed dust features in a representative mostly quies- 
cent sample of nearby (v < 3400 kms -1 ) early-type galax- 
ies using HST/WFPC2 observations. Unfortunately, this 
sample contains galaxies with a typical host magnitude 
~2 magnitudes fainter than in the UGC FR-I sample. It 
also has a fraction of lenticular galaxies which is about 
four times higher. We therefore cross-correlated the HST/- 
WFPC2 archive with the UGC catalogue to select all 
early-type UGC galaxies at v < 7000 kms -1 which have 
a similar absolute magnitude as the radio galaxies (i.e., 
M B < -20.4" 1 ), but which are not part of the UGC FR- 
I sample. The resulting sample of 52 galaxies is listed in 
Table and its general properties are shown in Figure [21 
Some of the non-FR-I galaxies host low-level activity in 
their nuclei, e.g., nuclear radio emission and optical line- 
emission. For clarity we therefore refer to this sample as 
the UGC non-FR-I sample instead of, for example, 
the quiescent sample. However, none of the UGC non- 
FR-I galaxies have jets on scales of tens of kiloparsec or 

1 Lyon-Meudon Extragalactic Database 

l |http:/ /leda.univ-lyonl ,fr| 



larger. The UGC non-FR-I sample turns out to match the 
UGC FR-I sample in velocity dispersion, recession velocity 
and morphologies (Figure [21 and Table |2J). 

The UGC FR-I and UGC non-FR-I samples also match 
in imaging depth and filter coverage. The entire UGC FR- 
I sample was observed in V- and /-band. Table [3 shows 
that 88% of the galaxies has V and/or /-band observa- 
tions, 54% has V- and /-band observations and 12% has 
only /?-band observations. Thus the filter coverage is sim- 
ilar for both samples. The many targets with both V- and 
/-band observations allow us to verify that the dust detec- 
tion rate and the derived dust properties do not depend 
on the filter used. For the UGC FR-I sample the expo- 
sure times for both V- and /-band images are typically 
460s (increasing up to 1400s in a small subset of the ob- 
servations). Table [3 shows that the exposure times for the 
UGC non-FR-I galaxies are similar. 

3. Dust Properties and Position Angle 
Measurements 

It is likely that every early-type galaxy contains dust due 
to mass loss from evolved stars. In addition, it can con- 
tain externally acquired dust from tidal interactions or 
accretion of galaxy satellites. Whether the dust obscura- 
tion is detected in optical observations depends on vari- 
ous factors, including sensitivity of the observations, dust 
morphology (e.g., smooth or clumpy) and orientation of 
dust with respect to the observer (e.g., face-on or edge-on). 
To ensure a homogeneous dust classification and measure- 
ment of dust properties for all samples we re-analysed the 
published HST imaging using the HST archive. We de- 
fine 'dust' as a localized depression in the stellar surface 
brightness. Thus we are insensitive to a global smooth 
dust component (e.g., following the stellar distribution) 
if it were present in early-type galaxies. Our approach is 
similar to that used in earlier HST imaging studies (e.g., 
van Dokkum & Franx 113351 Verdoes Kleijn et al. IT^l 
Tran et al. 15001) . 

The appearance of the dust structures in both ra- 
dio and non-radio galaxies varies from regular elliptical 
shapes, evoking the idea of inclined disks, to highly ir- 
regular structures, suggesting unsettled dust. We use the 
classification scheme for the UGC FR-I sample described 
in Verdoes Kleijn et al. (1999). It divides the morphol- 
ogy of the dust structure as projected on the plane of the 
sky into three bins. Class one, a 'dust lane', is a filamen- 
tary structure which passes through the centre and which 
is sufficiently regular to assign an orientation to it. Class 
two, a 'dust ellipse', is a dust structure with a circumfer- 
ence that resembles an ellipse. In some galaxies it is not 
possible to determine unambiguously if the dust has an 
ellipse- or a lane-like appearance. This is often due to the 
faintness and/or rather small angular size of the dust dis- 
tribution. We classify these galaxies as 'intermediate' (see 
Tables ^ and . Class three, 'irregular dust', is extended 
clumpy and/or filamentary dust which either is too ir- 
regular to assign an orientation or does not pass through 
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the nucleus and hence is not classified as a lane. Some 
galaxies which host a dust ellipse or lane also harbour a 
more extended, class three, irregular dust distribution (see 
Tables Hand EQ. We discuss in Sections and how the 
dust classification scheme, based on appearance, depends 
on viewing angle towards the dust distribution. 

We define the size of the dust as the largest linear ex- 
tent of the dust feature. This is the dust major axis for 
dust ellipses. We measure the ellipticity e of dust ellipses 
using their circumference. Dust classes one and two, the 
lanes and ellipses, are regular and extended enough to de- 
fine also a position angle PAd- For lanes it is the PA of 
the filament, taken as close as possible to the nucleus if 
any bending is present. For ellipses it is the PA of the dust 
ellipse major axis. For close-to-round ellipses (e < 0.1) we 
cannot determine a reliable PA of the major axis. The dust 
classification and dust properties for the FR-I and UGC 
non-FR-I samples are listed in Tables ^ an d El respec- 
tively. As a cross-check we compared our measurements 
of dust position angle and ellipticities with measurements 
available in the literature for various galaxies (Table EJ) . 
Ellipticities differ by ^ 0.05 and PA D differ by < 10°. 

We define the galaxy orientation as the position angle 
of the stellar isophotal major axis just outside the radius 
of the main dust distribution. For the UGC non-FR-I, 
the B2 and some 3CR galaxies, we measured the PA our- 
selves from isophotal fitting to the WFPC2 images, mask- 
ing dusty regions. The PA was taken from the literature for 
the remaining radio galaxies (see Tabled for references). 

For the radio galaxies we also require the position an- 
gle PAj of the jet axis close to the galaxy centre. For the 
UGC FR-I sample, we use the PAj as determined by Xu et 
al. (2000). These PAs are obtained from VLB A measure- 
ments, or, if those are not available, from the larger-scale 
VLA measurements. In cases where both VLA and VLBA 
measurements are available, the difference in PAj is 4° 
on average and always < 14°. For the other radio galaxies 
the PAj was taken from various studies (see Tabled for 
references). Given the difference in PAj between VLA 
and VLBA measurements and the published values for 
other jets, we estimate a typical error of 8° for all PAj. 

We are interested in the relative orientation of host, 
dust and jet if present. Thus we list in Tables Q and 0] the 
position angle difference APAdg between dust and galaxy 
axis and APAdj between dust and jet axis, both defined 
in the range [0°,90°]. 

4. Comparison of Dust Properties 

We first check whether the observed dust properties in our 
galaxy samples depend on distance. Figure |3 compares 
the radio galaxies (the UGC FR-I sample) and non-radio 
galaxies (the UGC non-FR-I sample). The top diagram 
shows that galaxies with and without dust detections are 
observed throughout the sampled volume. A similar range 
in dust sizes is detected at every distance except for small 
(< 100 pc) dust distributions which HST can resolve well 
only at small distances. We see no trend in morphology 



classification with distance. We conclude that the detec- 
tion of dust, the classification and its size do not correlate 
with galaxy distance for v < 7000 km s -1 . 

We detect dust in 48% ± 10% of the UGC non-FR-I 
galaxies. This detection rate is similar to the 40% ± 9% 
and 43% ± 8% detection rates from HST imaging for 
the samples of nearby early-type galaxies as compiled by 
van Dokkum & Franx l)I995|l — excluding galaxies with ex- 
tended radio structures in their sample — and by Tran et 
al. (2001), respectively. Both studies report an increased 
dust detection rate for galaxies with radio emission (com- 
pact and/or extended), but those rates are still less than 
the 90% ± 7% detection rate for our UGC FR-I sample 
of more powerful radio galaxies. At much lower radio lu- 
minosities (L ~ 10 19 ~ 21 W/Hz at 3.6cm), the radio lumi- 
nosity functions of nearby early-type galaxies with and 
without HST dust detection seem to become very similar 
(Krajnovic & Jaffe 2002). Hence, there is an increasing 
dust detection rate with increasing radio power. Many of 
the dusty non-radio galaxies have low-level active galactic 
nuclei, sometimes with radio cores (e.g., Ho et al. IT557I 
Tran et al. 2001). Thus the difference in detection rate is 
perhaps a trivial effect: it could reflect the simple fact 
that all radio galaxies have active galactic nuclei. The 
close connection between dust and nuclear activity then 
supports the idea that the central dusty ISM forms — not 
surprisingly — the fuel reservoir for the active nucleus. 

All three dust morphologies occur in both radio and 
non-radio galaxies. Dust features with a lane-like morphol- 
ogy appear in dusty radio and non-radio galaxies at simi- 
lar frequencies. In contrast, irregular dust distributions are 
more often seen in non-radio galaxies while radio galaxies 
harbour more often dust ellipses (see Table HJJ. Thus more 
powerful radio sources occur preferentially in galaxies with 
regular (perhaps more 'settled'?) dust. 

Lanes and ellipses span the same range in size. There is 
no clear size difference for these dust morphologies in radio 
and non-radio galaxies. The compact and/or faint features 
which could be either lanes or ellipses form the smallest 
structures at a given distance. The irregular dust distribu- 
tions tend to be the largest dust features. Larger-scale dust 
distributions tend to have a more irregular morphology. 

The middle panel of Figure |3| reveals no systematic 
difference between the host magnitude of each dust class 
for either the UGC FR-I or the UGC non-FR-I sample. 
Perhaps galaxies brighter than Mb ~ —21.5 do not con- 
tain dust with sizes below ~ 200pc, but this result is not 
statistically significant. 

The bottom panel of Figure shows a correlation be- 
tween dust orientation and morphology. The regular el- 
lipses align within ~ 10° with the galaxy major axis, but 
the filamentary dust lanes do not. This result holds for 
both radio and non-radio galaxies and is also seen in the 
sample of lower-luminosity early-type galaxies of Tran et 
al. (2001). A corollary of these results is that the mor- 
phological classification based on appearance cannot be 
due (only) to viewing angle, but reflects (at least partly) 
intrinsic differences between dust lanes and ellipses. 
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The dust lanes in radio galaxies have APAdg values 
over the complete range, i.e., APAdg ~ [0°,90°]. In con- 
trast, dust lanes in non-radio galaxies are within ~ 25° 
from either the minor- or major axis. It is not clear if 
this absence of large misalignment angles is simply due to 
the limited size of the sample. A similar absence was also 
seen by Tran et al. I|2001|l from WFPC2 imaging, but not 
by van Dokkum & Franx 1)1995(1 using WFPC1 imaging. 
At the same time (and regardless of morphology) smaller 
sized 1 kpc) dust structures tend to align within 30° of 
the galaxy major axis. Interestingly, only kpc-scale dust 
lanes, and the 120pc dust lane in NGC 2768, are roughly 
aligned with the galaxy minor axis. NGC 2768 has large- 
scale irregular dust with an extent of ~ 1.7 kpc roughly 
along the minor axis of the galaxy, i.e., parallel to its small 
dust lane. This irregular dust is associated with an ex- 
tended disk of emission-line gas which rotates around the 
major axis of the galaxy, i.e., orthogonally to the main 
motion of the stars (McDermid et al. 120041 Sarzi et al. 
2005). The similarity between radio and non-radio galax- 
ies in their relative orientations of dust lanes and ellipses, 
suggests, if taken at face value, that their orientations are 
not influenced by the presence of a radio jet. 

5. Relative Position Angles of Dust, Jet and Host 
in Radio Galaxies 

We now focus on the position angle of dust lanes and 
ellipses relative to (i) the radio jet axis and (ii) the galaxy 
major axis for the FR sample (see Table [Q. 

The left panel of Figure 01 is similar to the bottom 
panel of Figure |3 but now for the FR sample. It confirms 
the trends between the PA difference between dust and 
galaxy major axis APAdg and dust size of the subset of 
UGC FR-I galaxies in Figure |31 The segregation in the 
location of the data points as a function of apparent dust 
morphology is now less clear-cut. Dust ellipses still tend 
to align with the galaxy major axis (i.e., APAdg < 25°), 
but the alignment is not as tight as for the UGC FR-I 
sample by itself. Furthermore, there is one notable excep- 
tion at APAdg ~ 90° (3C 76.1, see Section for further 
discussion). Again, most lanes have large misalignments 
with the galaxy major axis (i.e., APAdg > 20°) and the 
largest lanes align roughly with the galaxy minor axis. 

The right panel of Figure 0] shows APAdg as a 
function of the PA difference between dust and jet axis 
APAdj- We note three special features. First, no radio 
jets are observed close to the dust longest axis, i.e., 
APAdj < 20°. Second, the data points are distributed 
roughly along a mirrored 'L' shape. Dust structures 
which are roughly aligned with the galaxy major axis 
(APAdg < 20°) have a wide distribution in relative 
angles with the radio jet (APAdj ~ [20° - 90°]). In con- 
trast, dust features which are misaligned from the galaxy 
major axis (APAdg > 20°) have a narrow distribution in 
PA differences with the radio jets (APAdj ~ [60° - 90°]). 
These dust structures, while misaligned with the galaxy, 
are in a very rough sense perpendicular to the radio 



jets. Third, dust ellipses and lanes have a different 
distribution of APAdj- Ellipses have a wide distribution, 
APAdj ~ 20° — 90°. In contrast, lanes have a narrow 
distribution: all lanes except one have APAdj ?=S 65°. 
The APAdj of the dust lane galaxy 3C 353, for which no 
APAdg is available (Tabled, is consistent with this. 

The FR sample contains 14 galaxies with dust fea- 
tures classified as 'ambiguous between ellipse and lane' 
(cf . Table ^ . APAdg an d APAdj can be measured for 
eight of them, and are plotted in Figure 0] Their distribu- 
tion among the lanes and ellipses supports the idea that 
in some cases these dust features are ellipses and in some 
cases they are lanes. 

What (selection) effect (s) could create the mirrored L 
shape in the right panel of Figure^f It is not unreasonable 
to assume that jets ploughing through a dusty medium 
would quickly destroy the dust. This could explain the 
absence of systems with APAdj < 20°. Furthermore, 
suppose that jets are roughly perpendicular to the dust 
plane intrinsically. If the dust structures are intrinsically 
circular and viewed close to edge-on, then we will observe 
e ~ 1 in combination with APAdj ~ 90° . If such a circu- 
lar disk is viewed close to face-on, we will observe e ~ in 
combination with a wide range of relative position angles, 
i.e., APAdj ~ [0°,90°]. The thin filamentary dust lanes 
might be interpreted as systems viewed relatively close 
to edge-on. Their high values of APAdj 60° are then 
qualitatively consistent with radio jets being intrinsically 
approximately perpendicular to the dust plane. 

6. Are the dust ellipses intrinsically circular? 

Before analysing the intrinsic orientations of dust ellipses 
and radio jets, we test if the observed distribution of e is 
consistent with circular or elliptical, thin or thick disks, 
assuming random viewing angles. This assumption is in- 
valid if (i) the disk detection rate or classification depends 
on disk ellipticity, or (ii) the radio selection criteria re- 
sult in a non-random distribution of jet viewing angles 
and there is a relation between disk and jet orientation. 
Scenario (i) is unlikely for the following reasons. We de- 
tect dust ellipses within almost the full range of possi- 
ble ellipticities (i.e., 0.01 to 0.86). There is no correlation 
between observed ellipticity and distance. There is a de- 
pendence of classification on distance for dust structures 
classified as 'intermediate' (i.e., between ellipse and lane). 
At distances D < 150 Mpc, ~ 14% of dusty galaxies in 
the radio and non-radio galaxy samples are classified as 
intermediate. These do not have a clear bias for e greater 
or smaller than 0.5. However, nine of the eleven galax- 
ies at D > 150 Mpc have an ambiguous dust morphol- 
ogy classification of which six have e ^ 0.5 (see Table 
At these large distances it becomes impossible with the 
current imaging to distinguish unambiguously between 
high-cllipticity ellipses and lanes in general. Thus we limit 
the analysis of disk shapes to galaxies at D < 150 Mpc. 
Scenario (ii) does not seem to be relevant for the UGC FR- 
I sample as it is selected on total 1400 MHz flux which 
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is dominated by the extended and unbeamed radio-lobe 
emission. The radio cores from VLA 1490 MHz measure- 
ments (FWHM ~ 1.5"-3.75") constitute always less than 
22% of the total 1400 MHz flux and typically ~ 6% (Xu 
et al. I2UU 0). Radio galaxies in the 3CR sample and B2 
sample are selected using total 178 MHz and 408 MHz 
fluxes, respectively fBennett Colla et al.. H970|l . At 

these lower frequencies the radio emission is expected to 
be even more lobe-dominated. 

We model dust ellipses as randomly-oriented, intrin- 
sically triaxial bodies with an intermediate-to-long axis 
ratio p and short-to-long axis ratio q < p < 1 . Thus an in- 
finitely thin circular disk has p — 1 and q — 0. We compare 
the observed distribution of e to the expected distribution 
for given p,q using the Kolmogorov-Smirnov (KS) test. 
Figure |5] shows the contour levels of equal KS probabil- 
ity in the plane of the axis ratios p and q, for the sample 
of 18 radio galaxies with dust ellipses at D < 150 Mpc. 
The maximum q allowed by the data is q = 0.14, given 
that the maximum observed ellipticity is e — 0.86. The 
contours are almost independent of q, suggesting that the 
data set is too small to put any further constraint on q. 
The KS probability is a function of p and maximizes for 
close to oblate (p ~ 0.74) disks. Similar results are ob- 
tained for the combined sample of 25 non-radio and radio 
galaxies at D < 150 Mpc with dust ellipses. In the analy- 
sis of the three dimensional configuration of jets and dust 
disks, we will consider three representative models: a thin 
circular disk (p = 1, q = 0), a thick circular disk (p = 1, 
q = 0.13) and a thin elliptic disk (p — 0.75, q = 0). 



7. Intrinsic orientation of dust ellipses and radio 
jets 

We now attempt to constrain the intrinsic orientation of 
the jet-axis relative to the dust systems classified as el- 
lipses in the FR sample. Figure illustrates the model 
parameters and coordinate system which we will use to 
describe the dust - jet system. The long, intermediate and 
short axes of the triaxial dust structure are chosen to lie 
along the X, Y and Z-axes, respectively. The radio-jet axis 
makes an angle #dj with the Z-axis and an azimuthal an- 
gle (f>Bj with the X-axis. The line of sight makes and angle 
#i os with the Z-axis (i.e., the disk inclination angle) and 
an azimuthal angle 4>\ os with the X-axis. Thus, the disk-jet 
model has in general six parameters. 



7.1. The relation between model parameters and 
observables 

The observations provide only two parameters of the 
dust-jet system as projected on the plane of the sky: the 
position angle difference between jet and disk APAdj 
and the disk ellipticity e. The ellipticity can be expressed 



in the model parameters as follows (Contopoulos 1956 
Franx lT§gg)> : 



{l-ef = 



with 



Vb' 



(1) 



a = (1 - q 2 ) cos 2 los + (1 - p 2 ) sin 2 <9 los sin 2 0i os + p 2 + q 2 , 
b= [(1 - q 2 ) cos 2 6>ios - (1 - p 2 ) sin 2 6> los sin 2 </> los - p 2 + q 2 } 
+4(1 - p 2 )(l - q 2 ) sin 2 9i os cos 2 6>i os sin 2 </> los . 

The position angle m i n of the minor axis of the dust 
structure relative to the projection of the Z-axis (i.e., 
short axis) can be written as: 



tan 20 n 



2T sin 4>i os cos ^i os cos 6» los 



sin 2 6»i os - T(cos 2 0i os - sin 2 <p\ os cos 2 0\ o 



;,(3) 



with triaxiality parameter T = (1— p 2 )/(l — q 2 ). The 
value of m in indicates the position angle of the minor 
(rather than the major) axis if: 



sign(tan 6 min ) = sign(sin0 los cos</>i os cos6>i os ). 



(4) 



The position angle 0j e t of the jet-axis relative to the 
position angle of the Z-axis can be expressed as (e.g., de 
Zeeuw & FranxUHSU): 



tan 



sin A</> tan #d j 



jet 



sin 9\ os — cos Acf) tan #d j cos 8i os 



(5) 



with Acj) = 4>\ os — 0D.I- The position angle difference 
0° < APAdj < 90° can then be expressed as: 



cos APA D j = sin(90° - APAdj) = |sin(e jot - K 



(6) 



Eq. shows that m in = 0° for a circular disk (p = 1 and 
q = 0), so that | sin6j C t| corresponds to sin(90° — APAdj). 
The expressions for e and APAdj then simplify to: 



e = 1 - |cos6>i os |, 
and 

tan(90°-APA DJ ) 



sin A</> tan dry j 



(7) 
.(8) 



sin 9\ os — cos Acf) tan 9u j cos 9i os 
Thus APA D j = 90° for any line of sight if 6> D j = 0°. 

7.2. Lower limits on #dj 

The angle between the line of sight and the radio jet is 
generally unknown. For an oblate dust structure (i.e., p = 
1 and O m in = 0°), the observed position angle difference 
with the disk merely confines the jet orientation to lie 
anywhere on the two 'jet-circles' which are the two great 
circles defined by the line-of-sight vector and the angle 
6j 0t (cf. Fig. Eland Eq.EJ there are two jet-circles because 
APAdj is an absolute value). The values of #dj and the 
difference in azimuthal angle A0 fix the jet location in 
the jet-circles. We will refer to #dj as the 'misalignment 
angle'. To observe simultaneously a given position angle 
difference APAdj < 90° and disk ellipticity e requires 
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a minimum misalignment angle Geometrically, this 

angle defines the latitude on a unit sphere which is tangent 
to the jet-circles. The jet axis intersects one of the two 
tangent points. The minimum misalignment angle can be 
expressed as (e.g., Schmitt et al. I2(J(J2|I 



sinflgj 1 = cos APA DJ sin6»i os . 



(9) 



Since most radio galaxies have APAdj 7^ 90° (Tabled, it 
follows that jets typically have 6>p] n > 0°, irrespective of 
our line of sight towards the system. The 9^ n for circu- 
lar/oblate dust disks are listed in Tabled The maximum 
flgj 1 required by the FR sample is 54° for 3C 449 and 
four of the 16 disks for which APAdj is available have 
6*Qj n > 40°. If the disks are circular or oblate, significant 
misalignments occur between jet axis and disk normal. 

7.3. The distribution of #dj-' thin disks 

We have lower limits on #dj for individual circular disk- 
jet systems. The next step is to constrain the full distri- 
bution of #d,t- This constitutes an inversion problem: we 
have to recover the distribution functions of misalignment 
angles 6*dj, azimuthal angles (f>uj, and line-of-sight direc- 
tions (#ios,<fros) from the sample of observational pairs 
(APA DJ ,e). We limit the analysis to the 18 dust ellipses 
at D < 150 Mpc which are consistent with a spheri- 
cally uniform random distribution of the line of sights (cf. 
Section El). We assume that 0dj has a uniform random 
distribution. The relative intrinsic orientation between jet 
and dust is characterized then by the misalignment angle 
distribution function Pdj(#dj). We interpret the distribu- 
tion function as a probability density function. Given the 
small number of observations, we explore three param- 
eterizations for Pdj instead of attempting a parameter- 
free recovery. They are illustrated in Figure [7] First, the 
observations are interpreted as being drawn from a step- 
function probability density distribution -Pdj(^dj) (angles 
in radians): 



< e m < e A , 

otherwise. 



(10) 



The angle 9 a is a free parameter. This step function, which 
we call model A, tests the hypothesis that the distribu- 
tion of #dj is peaked either at small or large misalignment 
angles or better agrees with a uniform distribution (i.e., 
B a = t/4). Quantitatively, a fraction 1 — 26 a/tt of the jets 
have misalignment angles #dj < 9 a, while the average 
#dj = 9 A- Thus a 9a = 7r/4 corresponds to a uniform dis- 
tribution in #d,t- Second, to test the hypothesis that the 
distribution of #dj does peak at intermediate misalign- 
ment angles (i.e., #dj = f/4) we analyse the data with 
the following two-step function, which we call model B: 



9b- Thus 9b = 7r/4 corresponds to a uniform distribution 
in #d,t- Third, we explore a Gaussian distribution for (9dj: 



fbj(0Dj)a0Dj = Cexp(-- - 



?DJ, 



< 6» D j < 



The Gaussian is truncated to the physically allowed region 
< 0dj < 7r /2 and normalized to 1 using constant C. 
Both mean /x and dispersion a (of the full Gaussian) are 
left free to vary in the ranges < fi < ir/2 and < a < it. 
This model, which we call model C, explores the presence 
of a peak not just at #dj = 7r/4 but anywhere and with a 
free width. 

First, we estimate the 9a that best fits the observa- 
tions using a Maximum Likelihood analysis. The goal is 
to maximize the likelihood L of observing our data set of 
APAdj and e: 



L = fj P(APA DJ)i) e i |^)d(APA DJ )de ) 



(13) 



i=l 



where the product is over the n Q b s observations. We ap- 
proximate P(APAdj, e\9 a) using a Monte-Carlo simula- 
tion. We draw a large number of realizations of random 
line of sights {9\ os , 4>\ os ), random azimuthal jet angles 0j O t 
and misalignment angles #dj from model A. Using Eqs J7J 
and (jHJ, these realizations are expressed in APAdj and 
e. A two-dimensional histogram of the realizations with 
n bin equal-sized bins approximates P(APAdj,j, u\9a) and 
hence the likelihood. The likelihood is maximized using 
an amoeba routine (Press et al. 11994(1 with 9a as the opti- 
mising parameter maximizing L. We choose the numerical 
simulation approach for two reasons. First, P(APAdj,c) 
is not available in analytical form in general. Second, the 
Monte-Carlo simulation allows us to incorporate easily ob- 
servational selection effects. For example, the three dust 
ellipses with e < 0.1 are too roundish to determine APAdj 
(see Section OJl . We account for this by selecting a subset 
of n mc pairs of (APAdj,c) from the Monte-Carlo realiza- 
tions which have e > 0.1. We performed simulations using 
n mc = 10 6 and 7ibi n = 10. The numerical error arising from 
this choice is negligible compared to the uncertainties in 
9a caused by the limited size and accuracy of the data set. 
Maximum Likelihood estimators are not guaranteed to be 
unbiased (e.g., Cowan 1998). We verified with simulations 
that the bias of the Maximum Likelihood estimator used 
here does not affect our results. 

The Maximum Likelihood analysis infers 9° A hs = 43° 
for the n bs = 15 observations with both APAdj and 
e > 0.1. We verify that the data are plausibly drawn 
from the best-fitting model by comparing the value of 
the Maximum Likelihood L™£* of the observations to 

r max 
sim 



PDj(0Dj)d0DJ = i 

Pdj(0dj)<Z0dj 



2 2~ 



L G?#DJ, 



of simulated data sets which were created using 
the best-fitting 9 A hs . These mock samples contain the 
ZL_£B<fl nl <ZL + £B same number of galaxies as the observations. We find 
^dfes* > ^sim m 60% of the mock samples and hence the 
hypothesis is statistically well accepted. 



4 2 
otherwise 



In this case a fraction of 1 — 29b/tt jets have 7r/4 — 9b/2 < 
9d.i < 7t/4+6'b/2 and the average #dj = 7r/4, regardless of 



We estimated the confidence levels on 9° A hs in a similar 
way. Figure |H1 (top-left) shows the cumulative distribution 
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function of 9 A m inferred from 100 mock samples which are 
created in three ways: 

1. In a Monte-Carlo fashion using the best-fitting 9 A hs in 
model A from the observations and the requirement 
that e > 0.1. 

2. 'Dithering' the observations within the measurement 
errors. 

3. By a bootstrap method, i.e., creating mock data sets 
of APAdj , e by drawing n oos independent data pairs 
(APArjj,e) from the n ODS observations. 

The first method samples the uncertainty of the best fit 
due to the finite number of observations under the as- 
sumption that the underlying model is correct. The second 
method samples the uncertainty due to the random obser- 
vational measurement errors. The third method samples 
a similar uncertainty as method 1 but with the observed 
data set itself used as an estimator of the underlying prob- 
ability distribution. 

All three methods infer median values for 6*^ m ~ 
40° — 45° in good agreement with the observations. The 
'dithering' method yields la, 2a confidence intervals which 
are significantly smaller than those of the other two meth- 
ods. This shows that the statistical spread in 6a due to the 
finite sample size is the dominant source of uncertainty for 
this model. The confidence intervals from the Monte-Carlo 
and bootstrapping method are more similar, although the 
Monte-Carlo method shows the largest 'probability wings' 
at small and large misalignment angles. Therefore we take 
the result from the Monte-Carlo method as a conservative 
estimate of the confidence levels. 

We conclude that the observed jet-disk relative po- 
sition angles and disk inclinations are consistent with a 
spherically random distribution of misalignment angles. 
The la upper and lower confidence levels indicate that 
the average #dj lies in the range 38° — 56°. 

7.4. The distribution o/~6>dj-' thick disks 

In Section we showed that the dust ellipses are not only 
consistent with thin circular disks (i.e., axis ratios p = 1 
and q = 0), but also with oblate (i.e., thickened) disks, 
with axis ratios p — 1 and q = 0.13. For such structures, 
Eq. Q simplifies to 

e = 1 - \J cos 2 9 los + q 2 sin 2 9 los (14) 

Thus, in comparison to thin disks, the thick disks appear 
rounder at any inclination 9\ os . As a result, the inferred in- 
clinations for the disks in the radio galaxy sample increase 
under the assumption of thick disks. Eq. @ then implies 
that the minimum misalignment angle 9^™ increases for 
each disk. We redid the Maximum Likelihood analysis as- 
suming dusty oblates with q = 0.13. The analysis for thick 
disks yields as expected a slightly larger, but very similar 
gobs _ 43 50 compared to thin disks. Also the confidence 
levels do not change significantly (see Figure ISJ. Thus the 
inferred distribution of jet misalignment angles does not 
depend critically on the assumed thickness of the disks. 



7.5. The distribution of9j)j: elliptic disks 

We showed in Section[|)]that the dust ellipses are also con- 
sistent with thin elliptic disks (i.e., p = 0.75 and q = 0) 
observed at random viewing angles. In contrast to circular 
disks, a APAdj < 90° can be observed for certain view- 
ing angles also for 9dj = 0°. However, only a fraction of 
the full range 0° < APAdj < 90° can be observed for 
certain e. Thus some dust ellipses still require a minimum 
misalignment angle 6>gj n > 0° in the case of an elliptic 
disk. We determined these angles numerically and they 
are listed in Table^next to those for the case of a circular 
disk. For three of the 16 dust ellipses 9^™ > 40° assuming 
p = 0.75. Thus, significant misalignments do occur also if 
the dust ellipses are thin elliptic disks intrinsically. 

We repeated the Maximum Likelihood analysis assum- 
ing that all dust ellipses are elliptic disks with p = 0.75 
and q = 0. The result is plotted in the middle- left dia- 
gram in Figure |SJ As expected, the best-fitting average 
misalignment angle 9 A hs = 40° is smaller than for thin cir- 
cular disks. The change falls within the la confidence level 
for circular disks. Thus, the best-fit to the distribution of 
jet misalignment angles does not depend critically on the 
assumed ellipticity of the disks within the accepted range 
p= [0.75,1]. 

7.6. The peak in the distribution of misalignment 
angles 

The absence of peaks at small or large misalignment an- 
gles resulting from using model A cannot rule out a best-fit 
distribution of misalignment angles which peaks at an in- 
termediate #dj ~ 45°. Thus it is worthwhile to explore 
model B which can test for such peaks. For circular disks, 
the analysis infers a best-fitting 9° B hs = 27°, i.e., a broad 
peak corresponding to^ 70% of the jets having misalign- 
ment angles 6*dj = [32°, 58°]. Figure |S] (top-right) shows 
the confidence levels using the three methods of parame- 
ter error estimation also used in Section IT"3l Although a 
broad peak is preferred, all methods indicate that neither 
a single #dj = 45° for all jets nor a distribution uniform 
in #dj (i-e., 9b = 45°) can be ruled out at more than the 
95% confidence level. 

Figure |H1 (middle-right) shows that the assumption of 
thick disks yields similar preferred 9b, while assuming thin 
elliptic disks yields 9b = 22°, i.e., a peak which is 5° 
smaller in width. However the upper 95% confidence level 
on 9b increases, while small 9b 10° are ruled out at 
larger confidence compared to thin circular disks. 

We explore model C, the truncated Gaussian, to test 
for the presence of a peak in the #dj distribution outside 
#dj = 45°. The best fitted truncated Gaussian from the 
Maximum Likelihood analysis indicates that at least half 
of the radio jets make an angle of ~ 35° or more with 
the symmetry axis of the dust disks (see Figure |j5J. The 
confidence intervals on the free parameters /1 and a are 
degenerate at large a. The reason is the #dj distribution 
asymptotically approaches a uniform distribution regard- 
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less of the value of /i. This precludes a direct interpretation 
of the confidence levels and hence we do not explore them. 

In summary, the jet misalignment angle distribution 
in galaxies with dust ellipses is consistent with having a 
peak around #dj ~ 45°. The width of this peak is not well 
constrained. The limited data set cannot rule out a spheri- 
cally random distribution of misalignment angles or a very 
narrow peak at more than 95% confidence. These conclu- 
sions do not depend critically on the assumed thickness or 
ellipticity. Regardless of the assumed model, typically at 
least half of the radio jets make an angle of 45° or more 
with the symmetry axis of the dust disks. 

7. 7. From axes to vectors 

In addition to the disk inclinations and dust-jet PA dif- 
ferences we can estimate the near side of dust disk and 
jet with respect to the observer in 11 of the radio galaxies 
(of which 10 are at D < 150 Mpc). This provides extra 
constraints on the distribution function of jet-disk mis- 
alignment angles. The flux and morphology asymmetries 
between the radio jets on both sides of the nucleus have 
been interpreted as being due to relativistic motion of 
the radio-emitting particles in intrinsically identical jets 
emerging on opposite sides from the nucleus (e.g., Laing 
et al. I1999[) . Thus the jet pointed nearest to the line of 
sight (i.e., the 'main jet') appears brighter as the radio 
flux is either Doppler boosted or less Doppler dimmed in 
comparison to the jet on the other side of the nucleus. 
The column density of stellar light obscured by the dust 
disk along the line of sight is larger for the near side of 
the disk than for the far side. Thus one expects a dif- 
ference in stellar surface brightness on both sides of the 
disk. Such an effect is indeed seen: in many cases, the 
major axis of the dust disk divides the disk into a more 
and a less obscured half. We identify the side with a lower 
stellar surface brightness with the near side of the disk. 
The surface brightness difference is not clearly present for 
disks with a low ellipticity. This is in qualitative agree- 
ment with the assumption that the disks are intrinsically 
close to circular in which case a lower ellipticity indicates 
a more nearly face-on disk for which the difference in the 
obscured amount of stellar light decreases. For 11 galaxies 
the observations reveal a clear radio flux asymmetry in 
the jets and a dust obscuration asymmetry which allows 
us to estimate the near side of disk and jet (cf. Tabled]). 

With the addition of near and far side information for 
jet and disk, the projections of the jet and dust axes in 
the plane of the sky become projections of a jet and dust 
vector. The main jet is projected against either the near 
side or far side of the dust disk (see Tabled . The absolute 
position angle difference between the two vectors can vary 
between 0° and 180° instead of the 0° < APA DJ < 90° 
for the dust and jet axis. We repeated the Maximum 
Likelihood analysis using this extra information for the 
ten galaxies at D < 150 Mpc. The 6 , ^ bs inferred for model 
A decreases by ^ 5°, indicating a slightly more peaked 



distribution of misalignment angles compared to the anal- 
ysis without near side information (see Figure HJ middle 
row). The inferred width of the peak at 6>dj = 45° for 
model B stays the same. 

For seven of the 11 galaxies more detailed modelling 
of the jet asymmetries at radio frequencies has been pub- 
lished (see Table^for references) . These studies report not 
only the near side of the jet but also an estimate of the 
viewing angle 6>jl to the main jet (cf. Table HJ. Under the 
assumption of circular disks, this viewing angle constrains 
the main jet to pass through either of the two points on 
each 'jet-circle' at which the circle defined by #jl inter- 
sects (cf. Section 17.20 . These four points correspond to 
two different #dj- For one point the main jet is projected 
against the near side of the dust disk while for the other 
point it is projected against the far side of the dust disk. 
Given our estimate of the near and far side of the dust 
disk, we can determine a unique 6*qj'° which is listed in 
Table m At least five of the seven misalignment angles are 
^radio -> 4qo^ connrrrun g that such large misalignments 
occur frequently. Special support for this comes from the 
galaxies with 10 which are much larger than the mini- 
mally required misalignment angle 6^\ n . In fact, the distri- 
bution of 0qj ° agrees especially well with the distribution 
inferred from models B and C (see Figure El . 

8. Intrinsic orientation of dust lanes and radio jets 

We need to define a three-dimensional orientation for dust 
lanes to constrain the jet misalignment angles in radio 
galaxies with dust lanes. The analysis in Section shows 
that the morphology of the dust ellipses is consistent with 
randomly oriented disks. The circumference of the dust 
lanes is too irregular to assign an ellipticity. It seems that 
the dust lanes are all viewed rather close to edge-on: the 
ratio of shortest to longest linear scale in the dust lanes is 
typically similar or smaller than the axis ratio of the most 
elongated dust ellipses in the sample. 

To obtain an upper limit on the average misalignment, 
we will assume that all dust lanes are exactly edge-on sys- 
tems. It follows from Eq. © that the minimally-required 
misalignment angle 6*^™ is maximal for 8\ os = 90°, i.e., an 
edge-on disk. Thus we expect to first approximation that 
the Maximum Likelihood analysis will infer an upper limit 
to the median of the distribution of #dj if ah lanes are as- 
sumed to be edge-on. We verified that this is indeed the 
case for the specific distribution of APAdj for the eight 
lanes at D < 150 Mpc. 

The upper limit to the median misalignment angle for 
dust lane radio galaxies is smaller than the median angle 
for dust disk radio galaxies (cf. Figure EJ. This conclusion 
is reached by comparing models A, B and C. The inferred 
angle 0b — 41° suggests that the data are not consis- 
tent with a peak in the misalignment angle distribution at 
#dj = 45°. On the contrary, the 9 a = 28° and the typical 
misalignment angle 6*dj ~ 20° from model C indicate that 
the observations are most consistent with such small mis- 
alignment angles. For model A the 9 a inferred for the dust 
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lanes falls outside the two-sided ~ 95% confidence region 
inferred for dust ellipses (cf . Figure |SJ) . The exact confi- 
dence level depends on the assumed thickness, ellipticity 
of the dust disks and whether additional assumptions on 
near side of the jets and disks are included. Similarly, for 
model B, the inferred 9b falls outside the two-sided ~ 68% 
confidence region for dust ellipses. 

The fact that dust lanes seem to be viewed close to 
edge-on prompts the question: where are the relatively 
face-on dust lanes? We explore several possibilities. It 
could be that close to face-on dust lanes are not classified 
by us as 'lanes'. At D < 150 Mpc there are five galax- 
ies with dust classified as intermediate between lane and 
disk. Three of these have axis ratios larger than 0.5 (i.e., 
roundish) and two have axis ratios less than 0.5 (i.e., flat- 
tened). Two of the three roundish galaxies have optical 
jets (3C 66B: Butcher et al. 115501 and 3C 78: Sparks et al. 
I1995fl . suggesting the jet is pointed close to the line of sight 
(e.g., Sparks et al. l2000[> . Thus at least these two galaxies 
are consistent with having dust structures which are close 
to face-on and having jets roughly perpendicular to them, 
i.e., face-on counterparts of dust-lane systems with per- 
pendicular jets. Furthermore, the two galaxies with axis 
ratios < 0.5 have APA D j = 68° and APA D j = 76°, i.e., 
the large angles typical for the dust lane galaxies. 

There are also five galaxies with dust classified as ir- 
regular. The dust in these galaxies typically shows several 
dusty filaments crossing the nucleus at several position an- 
gles. It is not implausible that these perturbed systems are 
the closer to face-on counterparts of perturbed dust lanes. 
The dust size distribution of the combined subset of galax- 
ies with either intermediate dust or irregular dust is simi- 
lar to that of dust lanes (see Figure0]and Table^l. In one 
galaxy with irregular dust (M87) an optical jet is seen. In 
fact, three of the four galaxies at D < 150Mpc with optical 
jets in our sample are found in galaxies with dust classified 
as intermediate or irregular. The one other optical jet is 
in a dust disk galaxy (NGC 3862). This galaxy has a very 
round dust structure (Table ^) . Thus one can speculate 
that this might also be a 'face-on dust lane' galaxy where 
the jet is perpendicular to the disk (see also Sparks et al. 
2000 ) . It is interesting to note that two other round dust 
ellipses (NGC 541 and UGC 7115) do not have optical 
jets. Thus these could then be the dust disk galaxies for 
which we inferred that the jets are not perpendicular to 
the disks. In conclusion there is evidence that many of the 
ten galaxies which contain intermediate or irregular dust 
could be the closer to face-on counterparts of the eight 
rather edge-on dust lane galaxies. 

9. Discussion 

We repeat the three questions posed in the Introduction 
which form the focus of our analysis. 

Are the properties of the fuel reservoirs different in 
radio galaxies and non-radio galaxies? 

In terms of dust morphology the answer is yes: non- 
radio galaxies more often display highly irregular dust dis- 



tributions. Perhaps this reflects the fact that the irregular 
fuel reservoirs are falling into the nucleus of the galaxy but 
are not feeding the black hole yet to cause nuclear activity 
and the formation of jets. In terms of reservoir orientation 
the answer is no. The orientations of dusty filaments and 
disks have similar distributions with respect to the optical 
isophotes for radio and non-radio galaxies. 

Is the dust orientation governed by the gravitational 
potential? 

Giant ellipticals are generally thought to be station- 
ary or at most slowly rotating triaxial structures. Such 
galaxies have two equilibrium planes for dust disks, i.e., 
planes where stable non- intersecting orbits exist (e.g., 
Merritt & de Zeeuw 1983 ). These planes are perpendic- 
ular to the short and long axis respectively. The angles 
between the minor axis and the short and long axis, 
respectively, depend on the line of sight and triaxiality 
T = (1 - p 2 )/(l - q 2 ) of the host galaxy (e.g., Franx et 
al. |l991). If we assume that ellipses and lanes are both 
intrinsically circular structures, then their minor axis and 
hence APAdg constrain the rotation axis assuming that 
the dust is settled. 

The top panel in Figure ITUI shows for the FR sample 
galaxies (at D < 150 Mpc) the observed cumulative dis- 
tribution functions of APAdg for dust ellipses and those 
predicted for randomly-oriented thin circular disks rotat- 
ing around the short and/or long axis as a function of 
host triaxiality. The bottom panel shows the correspond- 
ing probabilities from a Kolmogorov-Smirnov test that the 
observed distribution is drawn from a model distribution 
with a specific stellar host triaxiality. In computing the 
KS probabilities we added typical errors to the model dis- 
tributions. The errors have a skewed distribution because 
0° < APA D g < 90° by definition. Taking this into ac- 
count significantly influences the probability distribution. 
The APAdg of ellipses are only consistent with disks ro- 
tating around the short axis. The preferred host triaxiality 
T ~ 0.30 is close to oblate, in agreement with early studies 
of giant ellipticals (e.g., Franx et al. 1991). In this scenario 
of settled disks, one would expect to first approximation 
(i) dust ellipticities which are always larger than stellar 
host ellipticities and (ii) increasing host ellipticity for in- 
creasing disk ellipticity. This is exactly what is observed 

(Fig. CD. 

Examining the properties of individual galaxies instead 
of general sample properties , we note two radio galaxies 
which could be exceptions to the above scenario: NGC 
4261 and 3C 76.1. NGC 4261 is an unusual galaxy be- 
cause the stellar mean motion is around the major axis 
(Davies & Birkinshaw 1985J). This could indicate that the 
galaxy is prolate. In that case the dust disk would lie in a 
plane which is not perpendicular to the long-axis. Such a 
plane is dynamically unstable and one would expect that 
the gas would have settled down in another configuration. 
Thus it is more likely that the galaxy is triaxial but close 
to prolate. In that case the orientation of the dust disk 
and mean rotation axis of the stars are consistent with 
the picture in which the dust disk lies in a plane perpen- 
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dicular to the short axis of the galaxy, which is stable. 
The second stable plane in a triaxial galaxy, perpendicu- 
lar to the long axis, is more unlikely given that the disk 
is elongated in the direction of the galaxy major axis. 3C 
76.1 contains the only dust ellipse in the sample which 
is aligned with the minor axis of the galaxy instead of 
the major axis. Thus, this galaxy could be prolate with a 
settled dust disk in the stable plane perpendicular to the 
long-axis. Thus two out of 18 galaxies might be (nearly) 
prolate instead of (nearly) oblate. This fraction is consis- 
tent with the observed distribution of stellar rotation axes 
in large samples of ellipticals (Franx et al.[l991). 

To constrain the host shape for dust lane galaxies we 
start out with assuming that also dust lanes are thin cir- 
cular disks viewed at random viewing angles. Figure I1UI 
shows that the APAdg of dust lane galaxies is most con- 
sistent with equal numbers of dust lanes rotating around 
the short- and long axis in highly triaxial (T ~ 0.65) galax- 
ies. Such a high triaxiality appears inconsistent with other 
studies of the shapes of giant ellipticals (e.g., Franx et al. 
I1991JI . Furthermore, the similarity in large-scale axis ra- 
tios, optical luminosities and central stellar velocity dis- 
persion suggests that dust lanes and dust ellipses live in 
similar hosts (Figure 0). Figure ITUI shows that the data 
cannot rule out at high confidence that lanes occupy sys- 
tems with a typical triaxiality T ~ 0.2 — 0.3 but with lanes 
rotating around the long axis or around the short and long 
axis. However, as discussed in Section|Sl the assumption of 
random viewing angles appears implausible as dust lanes 
are more likely to be viewed relatively close to edge-on. 
Thus, we repeat the KS analysis but now we constrain 
the viewing angles to lie within 30° from edge-on. The re- 
sult is shown in Figure^] The hypothesis that all lanes 
are dust structures settled in the equilibrium plane per- 
pendicular to either the short or the long axis is ruled out 
with > 95% confidence. A model in which the dust lanes 
rotate in equal numbers around both axes cannot be ruled 
out at high confidence, but appears to be a poor fit to the 
observed APAdg a s well. 

The simultaneous differences between lanes and el- 
lipses in both their APAdg distribution and morphology 
support the alternative idea that the lanes are not simply 
settled, almost edge-on disks of dust. It is more likely that 
the dust lanes are perturbed disks. Models of a warped 
disk have been used to explain the filamentary dust mor- 
phology in NGC 4374 (M84; Quillen & Bower WMty and 
both the morphology and kinematics of the interstellar 
matter in NGC 5128 (Quillen et al. IUMjl . For NGC 4374 
and NGC 5128 the studies infer an angle between the 
line of sight and the symmetry axis of the warped disk 
of 68° and 60° respectively. In general, the perturbations 
and warps of the disks classified as 'lanes' could be non- 
transient if they are caused by triaxiality of the gravi- 
tational potential or constantly perturbing forces related 
to the active galactic nuclei and jets. The perturbations 
would be transient if the dust lanes are in the process 
of settling down to become flat circular disks (Steiman- 
Cameron et al. 1992), at which point we would classify 



them as dust ellipses. There is the following qualitative 
support for this. Some disks have in addition more irreg- 
ular dust at larger scales (e.g., NGC 315, NGC 4335 and 
NGC 4261). In particular the case of NGC 4335 is inter- 
esting as there the larger scale dust is attached to the disk 
and has an arc-like morphology very reminiscent of a set- 
tling disk (e.g., Verdoes Kleijn et al. 2002). Some of the 
most regular lanes have small APAdg (e.g., NGC 5127 
and NGC 5141). Furthermore, (some of) the five radio 
galaxies at D < 150Mpc with intermediate dust structures 
might be literally in between lane and ellipse. All have dust 
sizes < 550 pc which is typical for the settled dust ellipses, 
but their dust orientation is not aligned with the galaxy 
major axis. For the nine cases at D > 150 Mpc the clas- 
sification might be better described as 'ambiguous' due 
to large distance instead of intermediate (see Section EJ). 
Lastly, the hosts of dust ellipses seem close to oblate in 
general. In oblate galaxies the dynamical settling times are 
largest for lanes perpendicular to the equatorial plane. At 
the same time, the settling process is most likely dissipa- 
tional such that the dust structures become more compact 
as settling proceeds. Thus one expects larger lanes to have 
on average larger APAdg- This is observed (see Fig. |3J). 
Similarly, one expects large dust distributions to display 
closer alignment with the major axis at smaller radii be- 
cause of the decreased dynamical time scale. NGC 4335 
displays this behaviour, but it is not seen in general. 

It is reasonable to assume that the current dust dis- 
tribution was accreted by the galaxy before or around 
the time of the on-set of the AGN. The typical ages of 
FR I radio jets estimated from dynamical arguments is 
~ 10 8 yr and from radiative arguments ~ 10 7 yr (e.g., 
Eilek et al. I1997fl . Settling time-scales for the accreted 
gas are typically estimated to be a few orbital time-scales 
(e.g., Steiman-Cameron 1991). Orbital times are typically 
~ 2 x 10 6 (r/100pc) yr for giant ellipticals, where r is the 
radius for a circular orbit. (This estimate assumes that 
the orbital circular velocities are on the order of the cen- 
tral velocity dispersions as measured in the central few 
arcsec 2 .) The lifetimes of the AGNs seem a bit long for 
the dust still to be settling down. However, the uncer- 
tainties in both time-scales are too large to make firmer 
statements. 

Do jet/ AGN related forces affect the dust orientation? 

The dust/galaxy alignment dichotomy holds for giant 
elliptical galaxies in general, might they be radio or non- 
radio galaxies. A main result of our work is that for ra- 
dio galaxies the dust/galaxy alignment dichotomy is con- 
nected to a jet misalignment dichotomy in three dimen- 
sions. Filamentary dust structures, i.e., the dust lanes, 
tend to be perpendicular to jets and to have an arbi- 
trary orientation with respect to the galaxy major axis. 
Well-defined elliptical dust structures, i.e., the dust el- 
lipses, have an arbitrary orientation with respect to the 
radio jets and tend to be aligned with the galaxy major 
axis. We explore qualitatively two scenarios to explain the 
dust /galaxy /jet orientation dichotomy. 



12 



Verdoes Kleijn & de Zeeuw: A dichotomy in radio jet orientation 



Scenario I: the jets produce a torque in the hot X-ray 
gas. The torque is such that the X-ray gas tries to force 
the cooler dusty ISM in a plane perpendicular to the jet. 
Such a model was explored in detail for NGC 4374 by 
Quillen & Bower (1999). The stellar mass potential pro- 
duces a torque which tries to force the ISM in the stable 
settling planes of the galaxy potential. If the X-ray gas 
torque is stronger than the gravitational torque, the dust 
is forced in a plane perpendicular to the jet. This would 
produce the dust lanes galaxies. If the X-ray gas torque 
cannot overcome the gravitational torque the dust settles 
on non-intersecting orbits in the settling plane. This would 
produce the dust disk galaxies. It seems not implausible 
that an X-ray torque which is similar or stronger than the 
gravitational torque results in a perturbed disk with warps 
and filaments explaining the more irregular appearance of 
dust lanes compared to the dust ellipses (Quillen & Bower 
I1999fl . This scenario does not explain why the dust /galaxy 
misalignment is also seen in galaxies without radio jets. 
However, it is interesting to note that none of the dust 
lanes in non-radio galaxies have 20° APAdg < 80°. In 
particular, small kpc) dust lanes are oriented close to 
the galaxy major axis and all large (kJ kpc) dust lanes 
close to the minor axis. It might be that the small dust 
lanes in non-radio galaxies can be explained by dust struc- 
tures settled in planes perpendicular to the short axis in 
a mildly triaxial host. The large-scale dust lanes in non- 
radio galaxies might be still settling down towards this 
plane. In contrast to the non-radio galaxies, quite a few 
radio galaxies have misalignments at APAdg ~ 45° and 
quite a few of those have small sizes. This difference could 
indicate that dust lanes in radio galaxies are prevented 
from settling down by an additional force. 

Scenario II: dust is acquired externally (or from the 
outskirts of the galaxy) and the orientation of jet is aligned 
with the initial angular momentum of the dusty in- falling 
material. In this scenario the jet direction is determined 
during its formation by the angular momentum of in- 
falling material. A speculative example: it could be that 
an in-falling small galaxy brings in not only the dust and 
gas that we see, but also a central black hole. Initially the 
orientation of the dusty material then indicates the plane 
in which this black hole orbits the supermassive black hole 
of the giant elliptical. The angular momentum of the sub- 
sequent black hole binary or merger determines the orien- 
tation of the jets. Thus, initially the dust is expected to 
be perpendicular to the jet and to have a perturbed mor- 
phology because it has not settled yet. As time goes on 
the collisional dust and gas settles down in the galaxy po- 
tential thereby loosing its relation with the jet orientation. 
This scenario explains the dust/galaxy/jet dichotomy and 
the related differences in dust morphology. It also natu- 
rally explains why the dust/galaxy dichotomy is seen both 
in radio and non-radio galaxies. A weakness of the model 
is that the settling time-scales of the dust seem short com- 
pared to the age of the radio sources. 



10. Conclusions 

The main conclusions from the analysis of dust in 21 
nearby FR-I radio galaxies and a comparison sample of 
52 non-radio giant ellipticals are: 

— The detection rate of dust is higher in radio galaxies 
(~ 90%) than in non-radio galaxies (~ 50%). 

— The dust morphology can be classified as disk-like el- 
lipses, filamentary lanes and irregular. 

— Non-radio galaxies contain more often irregular dust 
distributions. 

— The dust morphology, size, orientation relative to the 
stellar host of ellipses and lanes is similar in radio and 
non-radio galaxies. 

— Filamentary dust, i.e., lanes, are generally misaligned 
with the galaxy major axis. Dust with a smooth ellip- 
tical appearance, i.e., dust ellipses, are roughly aligned 
with the galaxy major axis. 

— The dust lanes and ellipses live in similar hosts in terms 
of stellar luminosity, axis ratios and velocity disper- 
sions. 

The main conclusions from the analysis of the relative 
orientation (in two and three dimensions) of dust, host 
galaxy and radio jet in a sample of 47 radio galaxies are: 

— Dust ellipses are consistent with being relatively thin, 
circular or slightly elliptical disks rotating around the 
short axis of an oblate-triaxial gravitational potential 
with a triaxiality T ~ 0.3. 

— The dust lanes are most likely warped or per- 
turbed dust structures viewed rather close to edge-on. 
However, a scenario in which the dust lanes are settled 
structures in any of the two equilibrium planes of an 
oblate-triaxial galaxy cannot be ruled out completely. 

— The radio galaxies with dust with an apparent mor- 
phology classified as either irregular or intermediate 
between ellipse and lane most likely represent the close 
to face-on counterparts of the edge-on lanes. 

— The radio jets do not have a particular orientation with 
respect to the dust disks. The mean angle between 
the disk normal and jet (the 'misalignment angle') is 
inferred to be ~ 45°. 

— The radio jets tend to be perpendicular to the dust 
lanes, both in two and three dimensions. The upper 
limit to the mean misalignment angle is ~ 20° — 30°. 

— The above two results on jet orientations do not de- 
pend on the allowed range of ellipticity and thickness 
of the dust structures. 

— The use of independent information on the near side 
of jets and dust disks strengthen the above results on 
jet/dust orientation. 

There are two scenarios to explain the the 
dust/galaxy/jet orientation dichotomy for radio galaxies. 
Scenario I: the radio jets exert a torque on the nuclear 
dust in active galaxies which in some cases overcomes the 
gravitational torque causing the nuclear dust in a plane 
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perpendicular to the jets. Scenario II: the angular momen- 
tum vector of unsettled nuclear dust is initially aligned 
with the radio jet in active galaxies but this alignment is 
lost as dust settles in an equilibrium plane of the galaxy. 

To distinguish between these scenarios we deem the fol- 
lowing steps relevant. First, larger samples of dusty galax- 
ies, both active and inactive are required to confirm the 
dichotomy in jet/dust/galaxy orientation and its details. 
The current confidence levels on the existence of the di- 
chotomy are ^ 95%. Second, more detailed analysis of host 
galaxy properties, such as gravitational potential and X- 
ray gas pressure gradients, are required to obtain better 
estimates of the settling and life time of the dusty ISM. 
Third, combining such an analysis with detailed studies of 
the dusty and gaseous ISM, such as metallicity and two- 
dimensional kinematics, one might hope to constrain the 
origin of the dusty fuel reservoirs: stellar mass loss, cooling 
flows, or accretion of satellites. 
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Fig. 1. Central stellar velocity dispersion a versus total absolute blue magnitude Mb for early- type galaxies in the 
Uppsala General Catalogue at v < 7000 km s -1 . Circles indicate ellipticals, triangles are E/S0 transition objects and 
squares denote SOs. The filled symbols indicate the 18 radio galaxies of the UGC FR-I sample for which a a has been 
published. The horizontal and vertical solid lines mark the minimum a and galaxy luminosity observed for the radio 
galaxies. Magnitudes, dispersions and morphologies are taken from the LEDA catalogue. 
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Fig. 2. Central stellar velocity dispersion a, host absolute blue magnitude Mb, isophotal axis ratio (at a blue surface 
brightness of 25 mag/arcsec -2 ) and recession velocity V for the UGC non- FR-I sample of early-type galaxies (open 
symbols) and the UGC FR-I radio galaxies (filled symbols). Circles indicate ellipticals, triangles E/S0 are transition 
objects and squares denote SOs. 
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Fig. 3. From top to bottom: galaxy distance, absolute blue magnitude and position angle difference APAdg between 
host galaxy and central dust as a function of linear dust extent. UGC FR-I galaxies are represented by red symbols 
and UGC non-FR-I galaxies by black symbols. Circles indicate dust ellipses, triangles indicate dust lanes, squares 
indicate dust structures which have a morphology in between ellipse and lane and stars indicate irregular dust (see 
Sectional). The small dots indicate galaxies without dust detection plotted with a hducial linear dust extent of 20 pc. 
The typical relative error on dust size is 10%. The curve in the top panel indicates the minimum linear dust extent 
which can be resolved as a function of distance. In the bottom panel NGC 2768 is indicated because it has irregular 
dust with an extent of ~ 1.7 kpc roughly parallel to its small scale dust lane. See Section0]for further discussion. 
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Fig. 4. Left: PA difference APAdg between dust axis and host galaxy major axis as a function of dust size for the FR 
sample. Galaxies with a dust ellipse, a dust lane or an 'intermediate' morphology, i.e., in between lane and ellipse, are 
denoted by circles, triangles and squares, respectively. The subset of UGC FR-I galaxies have red symbols. The typical 
relative error on dust size is 10%. The trends in the full FR sample are similar to those observed in the UGC FR-I 
and UGC non-FR-I sample as shown in Fig. [31 Right: APAdg as a function of position angle difference between dust 
and radio jet APAdj. 
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Fig. 5. Contours of constant probability using a Kolmogorov-Smirnov (KS) test for the triaxial properties of the 
dust ellipses. The hypothesis is that the cllipticities e of the dust ellipses, observed in 18 F-I and FR-I/II galaxies at 
D < 150 Mpc, are caused by randomly-oriented triaxial bodies with an intermediate-to-long axis ratio p and short-to- 
long axis ratio q. Each contour is labelled with its KS probability. The distribution of e in the radio galaxies requires 
relatively thin (q < 0.14) disks and is most consistent with close to circular disks (p >^ 0.7). A similar result is 
obtained for the combined sample of 25 non-radio and radio galaxies with dust ellipses at D < 150 Mpc. See Sectional 
for details. 
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Fig. 6. The coordinate system used to describe the dust - jet system. The dust structure is assumed to be centred 
on the galaxy nucleus and to have a triaxial shape. The long, intermediate and short axis of the dust distribution lie 
along the X, Y and Z-axis respectively. The jet axis is indicated by the dotted arrow and makes an angle #dj with 
the Z-axis and an azimuthal angle 0dj with the X-axis. The line-of-sight direction, indicated by the dashed arrow, 
makes and angle 9\ os with the Z- axis and an azimuthal angle 4>\ os with the X-axis. The spherical angle labelled 0j ot 
indicates the angle between the jet and the short axis projected in the plane of the sky. See Section mi for formulae 
relating these quantities. 




Fig. 7. The three probability distributions of the jet / dust disk misalignment angle 6*dj assumed to underlie the 
observations. Left: model A, the single-step function as described by Equation^] The width and height of the level 
at low 8dj depend on the free parameter 9 a- The integrated probability underneath this level is 1 — 29 a /tt- This 
model is used to test if the observations imply a peak in the #dj distribution at low or high 6*dj- Middle: model B, 
the double-step function as described by Equation ^2 The width and height of the central level are set by 9b and is 
centred on #dj = 7r/4. The integrated probability under the central level is 1 — 29b/tt- Model B tests if the observations 
are consistent with a central peak or dip in the distribution of #dj- Right: model C, a truncated Gaussian normalized 
to 1 (see Ea. I12|l . Both the mean < /i < ir/2 and dispersion < a < ir are free parameters. This model tests for the 
presence of a peak anywhere in the range < 9uj < tt/2. See Section 1731 for details. 
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Fig. 8. Confidence levels on the best-fitted parameter 8a for model A (left column)and 8b for model B (right column) 
of the jet-dust misalignment angle, #dj, distribution. Plotted are cumulative distributions of the best- fitted free 
parameter, 8a and 8b, for various models of errors and dust morphology. The first four plots are models for the dust 
ellipses, the bottom two for the dust lanes. Top: thin circular disks are assumed (p = 1 and q = 0). To estimate the 
error on 8a the bootstrap (solid), Monte-Carlo (dashed) and dither (dotted) method are used (see Section EJ. The 
two-sided 68.27% and 95.45% confidence levels around the median are indicated by the horizontal dotted lines. The 
vertical lines indicate the best-fitted 8a and 8b for the observations. Middle: similar to the top model, but now for 
oblate 'thick' disks (p = 1, q = 0.13, dashed curve) and thin elliptic disks (p = 0.75, q = 0, dotted curve) and for thin 
circular disks (p = 1, q = 0, solid curve). The latter model uses additional information on the near side of jet and 
disk which is available for 10 of the 15 galaxies. All curves use the bootstrap method for error estimation. Bottom: 
similar to top diagrams, but now for eight lane galaxies. All models assume edge-on disks. The three methods of error 
estimation are indicated by the same line style as in the top diagrams. The main conclusion of the six diagrams is that 
dust lanes are most consistent with smaller misalignment angles than dust ellipses. The difference in misalignment 
angle is more clearly inferred by model C as shown in Figure See Sections [7] and [S] for a detailed discussion. 
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Fig. 9. Cumulative distributions of the jet-dust misalignment angle 6*dj for the radio galaxies with dust ellipses (top) 
and dust lanes (bottom). Top: the three thin curves indicate the cumulative distributions obtained from the Maximum 
Likelihood analyses of the observed distribution of dust ellipticities e and dust-jet position angles differences APA dj 
from 15 FR-I and FR-I/II radio galaxies (see Section UJ. The solid curve corresponds to model A (i.e., single-step 
function, cf. Equation HJ^i . the dotted curve corresponds to model B (i.e., two-step function, cf. Equation ITT|l and the 
dashed curve for model C (i.e., the truncated Gaussian distribution, cf. Eouation ll2|l . The models assume thin circular 
dust disks and take into account the information on the near side of disk and jet, which is available for 10 radio 
galaxies (see Section IT"T)l . The thick solid curve indicates the cumulative distribution for the seven radio galaxies for 
which estimates of individual $dj are available (see Section ITT)! . The shape of the latter curve does not depend much 
on the choice of (?dj for 3C 296 within its allowed range 30° < #d,t < 73° (we used the lower limit). Bottom: same as 
top panel, but now for the set of eight dust lanes. The lanes are assumed to be systems which are viewed edge-on (see 
Section[HJ). Correspondence between line style and model is similar to that in the top plot. The thick curve, present in 
the top plot, cannot be determined for dust lanes as e and individual estimates of 6*dj are not available for lanes. The 
main conclusion is that jets of radio galaxies with dust lanes have on average smaller misalignment angles than those 
in radio galaxies with dust ellipses. See Section |H1 for details. 
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host triaxiality T 



Fig. 10. Top: the cumulative distribution of APAdg- The two thick lines show the observations for the FR sample 
galaxies at D < 150 Mpc for dust ellipses (solid line) and dust lanes (dotted line). The thin dotted lines show the 
expected cumulative distribution of APAdg for circular disks rotating around the short axis (top series) and around 
the long axis (bottom series) assuming random viewing angles. The thin solid lines are the combination of these two 
distributions. They indicate the cumulative distribution for an equal number of disks rotating around the short and 
long axis. The nine curves in each series indicate the expectation for a stellar host triaxiality of 0.1 (top) increasing 
to 0.9 (bottom) in steps of 0.1. Bottom: the Kolmogorov-Smirnov (KS) probabilities that the observed APAdg 
distributions in the top plot are drawn from modelled distributions as a function of the stellar host triaxiality T. 
The computation assumes that the dust is a thin circular disk at random viewing angles and takes into account 
observational errors. The KS probability is significant at certain T only for dust ellipses rotating around the short axis 
(solid line), lanes rotating around the long axis (dotted line) and lanes rotating in equal numbers around the short 
and long axis (dashed line). Figure IT21 shows a similar analysis for dust lanes assuming viewing angles which are closer 
to edge-on. See Section |3] for further details. 
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ellipticity dust 



Fig. 11. The ellipticity of the stellar host as a function of dust ellipticity for the dust disks in the UGC FR-I sample. 
The host ellipticity is taken from HST imaging at a radius just outside the main dust distribution, using the results 
in Verdoes Kleijn et al. I|199 9). As expected for settled dust rotating around the short axis, the dust ellipses have 
ellipticities which are (i) systematically larger than host ellipticities and (ii) increasing for increasing host ellipticity. 
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Fig. 12. Similar plots as Figure but now for a model in which the thin circular disks are viewed within 30° 
from edge-on instead of fully random viewing angles. This model is more plausible for dust lanes. The Kolmogorov- 
Smirnov probabilities for lanes are only significant for lanes rotating around the short and long axis of the triaxial 
host. Nevertheless, this model does not provide a good fit to the data for any triaxiality T. In other words, it appears 
unlikely that the dust lanes are consistent with being settled structures in any of the two equilibrium planes in a 
triaxial host, but the data set is too small to rule out such a scenario completely. 
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Table 1. Dust properties for the FR sample of radio galaxies. Properties of dust detected in the FR sample, which includes 
the 21 UGC FR-I radio galaxies as a subset (see Section |2j. The latter are listed above the double horizontal lines. No dust 
is detected in two UGC FR-I galaxies: NGC 741 and NGC 2892. Col.(l): target name. A * indicates radio sources which are 
sometimes classified as FR-Is (e.g., de Koff et al. '20001 and sometimes as FR-IIs (e.g., Chiaberge et al. l20U2t or as intermediate, 
i.e., FR-I/II. Col. (2): distance to the galaxy taken from Faber et al. 115551 or LEDA, except for NGC 5128 (Israel 115551 . Col. (3): 
dust morphology classification (see Section |3] for their definition). 1= dust lane, 2— dust ellipse, 3= irregular dust. A 1-2 
classification indicates an intermediate morphology classification: they could be either lanes or ellipses. Col. (4): longest linear 
extent of the dust which has a typical relative error of 10%. Entries marked with a have an additional irregular extended 
component. Col.(5): ellipticity for dust ellipses. The error is given in between the brackets. We measured the ellipticity e from 
the HST imaging except for 3C402N, which was taken from Schmitt et al. (2002). Col. (6): dust position angle with its error in 
between brackets. This is the PA of the major axis for dust ellipses and of the longest dust extent for dust lanes. An 1 indicates 
that the dust feature was too faint to measure its position angle reliably. Col. (7): the position angle difference between the 
dust and galaxy major axis. The error is listed in between brackets. A 2 indicates that the galaxy is too round to determine its 
position angle. Col. (8): position angle difference between jet axis and dust ellipse major axis or longest axis of the dust lane 
with the error in between brackets. Col. (9): minimum misalignment angles 8^) n between the jet axis and the normal of the 
dust disk. The first and second value are obtained assuming thin (q — 0) circular disks (p = 1) and elliptic disks (p = 0.75), 
respectively. Col. (10): side of dust disk against which the main jet is projected; near/far indicates the side nearest or farthest 
from the observer, respectively. See Section l7.7l for details on its determination. Col. (11): viewing angle to the main jet (i.e., the 
angle between the line of sight and the jet) as inferred from jet asymmetries at radio frequencies. This angle is taken from the 
literature (see references). Col. (12): The jet misalignment angle as inferred from the viewing angle to the main jet. Col. (13): 
reference for PA of jet and galaxy. 1: de Koff et al. (2000); 2: Martel et al. <|1999fl : 3: Capetti et al. I|2()()()^ 4: Schmitt et al. 
ijUjJJJ; 5: de Koff et al. iflTMjl : 6: de Juan et al. (fT55H)l : 7: Guthrie (fT9"8"0|l (we measured the jet PA from Fig. 1); 8: Dufour 
et al. CnZnj; 9: Schreier et al. iflTMJl : 10: VIMOS imaging (unpublished); 11: Verdoes Kleijn et al. JT5J5J; 12: Verdoes Kleijn 
et al. (1277021 : 13: Leahy & Perley IT55TY References for 8 JL : 13: Laing et al. (12110111 : 14: Laing et al. JT555J; 15: Feretti et al. 
JT555J; 16: Hardcastle et al. (tT59Tll 



LcLI 


1 ' 


TT~ 

niorpn 


size 


6 


rAD 


iirADG 


iirADJ 


/jmin 
U D.J 


T-\ 

S1QG 


f JL 


/iradio 


j 

ret 




(Mac) 




(pc) 




(°) 


(°) 


(°) 


(°) 
\ ) 




(°) 


(°) 




NGC 193 


58 


i 


840 1 




0(6) 


58(7) 


77(10) 










ii 


NGC 315 


68 


2 


830 1 


0.77(0.02) 


40(2) 


1(4) 


90(8) 


0/0 


far 


37 


40 


11,13 


NGC 383 


65 


2 


2350 


0.23(0.03) 


138(2) 


6(3) 


24(8) 


36/0 


far 


52 


47 


11,13 


NGC 541 


73 


2 


640 


0.09(0.05) 
















11 


UGC 1841 


85 


1-2 


330 


< 0.5 
















11 


NGC 2329 


77 


2 


750 


0.32(0.05) 


174(3) 


3(4) 


24(9) 


42/29 


near 






11 


NGC 3801 


43 


1 


4430 1 




24(2) 


83(4) 


89(8) 










11 


NGC 3862 


84 


2 


610 


0.01(0.02) 
















11 


UGC 7115 


91 


2 


570 


0.05(0.03) 
















12 


NGC 4261 


30 


2 


240 1 


0.54(0.03) 


163(1) 


8(5) 


76(8) 


12/4 


far 






11 


NGC 4335 


62 


2 


750 1 


0.59(0.10) 


158(2) 


2(5) 


79(8) 


10/3 








11 


NGC 4374 


15 


1 


990 




79(10) 


49(11) 


80(13) 










11 


NGC 4486 


15 


3 


560 


















11 


NGC 5127 


64 


1 


950 




48(2) 


20(5) 


70(8) 










11 


NGC 5141 


71 


1 


550 




88(8) 


23(11) 


76(11) 










11 


NGC 5490 


70 


1-2 


170 


k, 0.5 


143(4) 


38(6) 


68(9) 










11 


NGC 7052 


55 


2 


1080 


0.70(0.02) 


65(1) 


1(1) 


42(8) 


45/40 


far 


79 


47 


11,14 


UGC 12064 


68 


2 


1200 


0.54(0.05) 


166(5) 


6(5) 


25(9) 


54/43 


near 


83 


72 


12,15 


NGC 7626 


47 


1 


230 




175(4) 


4(6) 


57(9) 










11 


NGC 5128 


3.4 


1 


9200 


> 0.5 


120(3) 

i 


85(4) 


67(9) 










8,9,10 


3C15* 


292 


1-2 


985 


< 0.5 














2 


3C29 


180 


1-2 


486 


< 0.5 


i 














2 


3C40* 


72 


2 


146 


0.53(0.07) 


77(7) 


25(8) 


62(11) 


24/17 


far 






2 


3C75 


96 


1-2 


194 


> 0.5 


36(5) 


65(5) 


76(9) 










2 


3C76.1 


128 


2 


892 


0.65(0.05) 


39(4) 
i 


89(5) 


72(9) 


17/11 








13 


3C78 


116 


1-2 


469 


< 0.5 














2 


3C83.1 


102 


2 


1020 


0.86(0.03) 


168(2) 


2(4) 


82(8) 


8/6 








1,2 


3C84 


72 


3 


12060 1 


















2 


3C88* 


121 


1-2 


550 


< 0.5 


l 














1,2 



28 Verdoes Kleijn & de Zeeuw: A dichotomy in radio jet orientation 



Table 1. continued. 



target 


D 


morph 


size 


e 


PA 


APA 


APA 




side 


n 

fjL 


/iradio 


ref 




{wipe) 








( ) 


( ) 


V ) 


( ) 




\ ) 


\ ) 






yo 



z 


OOO 


n 71 fn n^ 

U. f 1 (U.UO J 


1 H7fd^ 

10 C(4J 

1 


1 01 A\ 
1Z(4J 


ou(y ) 


00, lOO. 
zy / zo 


— p 


~ OO 


•?n 7i 

ou — ( 


1 1 R 




A QQ 

4oo 


1 o 
1-z 


148U 


<; U.o 














1 X. 

1,5 


opqi 7 


1 /in 

l^U 


q 
o 


1 1 1 7fl 
111 f U 


















o 
z 


oUooo 


iiy 


Q 

o 


I 7^n 

I I ou 


















1 9 

1 ,z 


oUoOo 


1 00 
Izz 


1 
1 


140U 




1 /I 8^7*1 
148^ / j 


2 


04(11 J 










1 9 


Qf^ 1 A HO AT 


1 m 
1U1 


z 


1 flQfl 

1U8U 


n rq/'h nrA 
U.Oo^U.UO j 


oo(o ) 


3(5) 


RQ f Ci\ 

0o(y ) 


Or; /m 
z0/ ly 








A 

4,0 




1 HQ 
lUo 


Q 

o 


COR 
ozD 


















o 
z 




30/1 
OZ4 


1 
1-Z 


msn 

oUOU 


u.o 


U(0) 




/ y (yj 










o 

z 


A HQ /I Q 

4U-Uo.4o 


ZU 1 


1 

1-z 


40U 


fC U.O 


80(0 ) 


08(0 J 


0y(y ) 










4 


do nn^/i-LO^ 

JjZ UUo4-|-ZD 


1 OS 
IZo 


o 
z 


1 01 n 

1Z1U 


U. / Z^U.UO J 


1 RV\( A\ 
10U(4 ) 


(\(A\ 
U(4J 


o / (yj 


00 /1 7 
ZZ / 1 i 


lar 


PI A 
04 


OA 
Z^± 


0,14 


DO HQ1 QO 


O/l 8 


o 


1 ^/in 

104U 


U.Z8(U.U4 J 


1 oof 1 nA 

IZZ^IU J 


oo<4 1 
zzyi i j 


88^1 *?\ 


1 /n 
l/U 


near 


an 

oU 


OO 


0,14 


P.9 1 9^^J_9S 
r>Z IZOD-rZo 


on 
yu 


z 


oOU 


n ai (c\ i n\ 

U.81 (U.1U J 




/I ('/I A 


45(8) 


A P. 1 AO 
40 / 4o 








^ 7 






1-2 


580 


^05 


on 0"l 


1 3(4 s ! 


30f1 3^ 










o 


B2 1346+26 


253 


1 


6710 1 




134(7) 


62(11) 


71(11) 










3 


B2 1357+28 


251 


1-2 


520 


^ 0.5 


95(5) 


2(10) 


85(9) 










3 


B2 1457+29 


588 


1-2 


5000 


£ 0.5 


36(4) 


78(5) 


61(9) 










3 


B2 1525+29 


261 


1-2 


540 


;> o.5 


148(5) 


28(5) 


62(9) 










3 


B2 2335+26 


120 


2 


1160 


0.25(0.05) 


6(9) 


24(9) 


61(12) 


19/0 


far 






3 



Verdoes Kleijn & de Zeeuw: A dichotomy in radio jet orientation 



29 



Table 2. Host morphologies. Relative distribution over host morphology (taken from the LEDA catalogue) for the 
UGC, UGC FR-I, UGC non-FR-I samples of galaxies. Col. (2): total number of galaxies in the sample. Col.(3)-(5): the 
fraction of galaxies within each sample for the three morphology classes and the Poissonian error. 



Sample 


N 


E 


E/SO 


SO 


comments 


UGC total 
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0.34 ±0.02 


0.16 ±0.01 


0.50 ±0.02 


v < 7000 kms" 1 




389 


0.49 ± 0.04 


0.17 ±0.02 


0.34 ±0.03 


v < 7000 kms" 1 and M B < -20.4 


UGC FR-I 


21 


0.71 ±0.18 


0.19 ±0.10 


0.10 ±0.07 




UGC non-FR-I 


52 


0.71 ±0.11 


0.12 ±0.05 


0.17 ±0.06 





Table 3. UGC non-FR-I sample. General properties of the 52 UGC non-FR-I galaxy sample which forms the compari- 
son sample for the UGC FR-I sample of radio galaxies. Galaxy data are taken from the LEDA catalogue. Cols. (l)-(3): 
galaxy name, Hubble type and absolute blue magnitude. Cols.(4)-(5): WFPC2 filter name and exposure time. Col. (6): 
HST program number. 
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6357 


5846 


E 


-21.2 


F555W 


2200 


5920 


3842 


E 


-22.2 


F555W 


1600 


6587 








F814W 


2300 


5920 


3894 


E-S0 


-21.0 


F547M 


503 


5924 


5982 


E 


-21.5 


F555W 


1000 


5454 


4073 


E 


-22.4 


F555W 


1600 


6587 








F814W 


460 


5454 


4125 


E 


-21.3 


F555W 


1000 


6587 


6211 


SO 


-20.9 


F606W 


500 


5479 








F814W 


2000 


6587 


6703 


E-S0 


-20.8 


F814W 


320 


5999 


4168 


E 


-20.7 


F547M 


460 


6837 


7318A 


E 


-20.9 


F569W 


1600 


6596 








F702W 


1000 


6357 








F814W 


1000 


6596 


4365 


E 


-20.9 


F555W 


1000 


5454 


7562 


E 


-21.1 


F555W 


2200 


6554 








F814W 


460 


5454 








F814W 


2200 


6554 


4406 


E 


-20.8 


F555W 


1000 


5454 


7619 


E 


-21.8 


F555W 


2200 


6554 








F814W 


460 


5454 








F814W 


2200 


6554 


4472 


E 


-21.4 


F555W 


460 


6673 


7785 


E 


-21.3 


F555W 


800 


6587 








F814W 


460 


6673 














4473 


E 


-21.6 


F555W 


1800 


6099 




















F814W 


2000 


6099 
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Table 4. Dust properties UGC non-FR-I sample. Properties of the dust detected in UGC non-FR-I sample galaxies. 
Col. (2): galaxy distance from Faber et al. (1989) or LEDA, except for NGC 4526, which is taken from Tonry et 
al. (|2001). Col. (3): dust morphology (see Section [3] for definitions). 1= dust lane, 2= dust ellipse, l-2=either lane 
or ellipse and 3=irregular dust. Col. (4): longest linear extent of the dust which has a typical relative error of 10%. 
Entries marked with a 1 have an additional irregular extended component. Col. (5): ellipticity of dust ellipses, with 
the measurement error in between brackets. Col. (6): dust position angle and its error in between brackets. This is the 
PA of the major axis for dust ellipses and of the longest dust extent for dust lanes. Col. (7): position angle difference 
between the galaxy major axis and dust longest axis and its error in between brackets (see Section |3J). 



NGC 


D 

(Mpc) 


morph 


size 

(PC) 


e 


PA D 

(°) 


APAdg 

n 


910 


69.23 


1-2 


160 


k, 0.5 






1129 


69.67 


2 


670 


0.85(0.10) 


0(2) 


0(3) 


1161 


25.80 


2 


2960 


0.53(0.05) 


135(5) 


2(7) 


1497 


81.69 


3 


5650 








2258 


53.12 


3 


470 








2768 


20.43 


1 


120 1 




53(8) 


84(9) 


2872 


41.65 


2 


150 


0.46(0.06) 


171(5) 


2(5) 


2911 


42.40 


1 


1080 1 




63(5) 


85(11) 


UGC3426 


53.73 


3 


3040 1 








3516 


34.95 


3 


2290 








3894 


42.95 


1 


1490 1 




109(4) 


14(4) 


4125 


26.48 


1 


170 1 




108(10) 


13(18) 


4406 


17.77 


1-2 


50 


^ 0.5 






4472 


17.77 


1 


210 




136(10) 


26(11) 


4494 


9.267 


2 


70 


0.50(0.05) 


23(6) 


1(6) 


4526 


16.90 


2 


2350 


0.80(0.05) 


127(5) 


5(7) 


4552 


17.77 


1-2 


40 1 


& 0.5 






4589 


40.40 


1 


2310 1 




171(10) 


89(11) 


4952 


79.11 


2 


190 


0.54(0.03) 


30(5) 


2(5) 


5252 


89.75 


3 


2510 








5322 


22.15 


2 


370 


0.87(0.05) 


88(5) 


2(5) 


5813 


31.15 


3 


1690 








5846 


31.15 


3 


2840 








7318A 


88.80 


3 


3790 








7785 


60.05 


3 


1610 
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Table 5. Comparison between published and our measurements of dust ellipticities and PAs. Col.(2)-(3): ellipticity 
as determined by us with the error in between brackets and as reported in the literature, respectively. Col.(4)-(5): 
position angle of the dust as determined by us (with error again between the brackets) and as reported in the literature, 
respectively. Col.(8): references 1: de Koff et al. PTO : 2: Capetti & Celotti (fT^jl ; 3: Capetti et al. (j2UUU|> : 4: Schmitt 
et al. 1(2003); 5: van der Marel & van den Bosch l(TUMjl : 6: de Ruiter et al. <|2UU2|l : 7: Dufour et al. l(T§75|) . The reason 
for the significant difference in the PAd for B2 1346+26, B2 1525+29 as measured by de Ruiter et al. l|2002j) and us on 
the one hand and by Capetti et al. |2000) on the other hand is unknown, but the de Ruiter et al. (|2002p measurements 
are deemed more accurate (Capetti private comm.). 



name 


c 


£rcf 


PA D 

n 


PA rcf 

n 


ref 


NGC 383 


0.23(0.03) 


0.18 


138(2) 


135 


2 


NGC 4261 


0.54(0.03) 


0.58 


163(1) 


165 


2 


NGC 7052 


0.70(0.02) 


0.65 


65(1) 


65 


5 


3C449 


0.54(0.05) 


0.50 


166(5) 


169 


1 


NGC 5128 






120(3) 


122 


7 


3C83.1 


0.86(0.03) 


0.91 


168(2) 


171 


1 


3C296 


0.71(0.05) 


0.71 


157(4) 


160 


1 


3C353 






148(7) 


166 


1 


3C402N 






58(3) 


55 


4 


4C-03.43 






86(5) 


73 


4 


B2 0034+25 


0.72(0.06) 




160(4) 


160/160 


3/6 


B2 0915+32 


0.28(0.04) 




122(10) 


125 


6 


B2 1256+28 


0.81(0.10) 




178(3) 





6 


B2 1339+26 






0(10) 


0/175 


3/6 


B2 1346+26 






134(7) 


0/134 


3/6 


B2 1357+28 






95(5) 


95/77 


3/6 


B2 1457+29 






36(4) 


30/50 


3/6 


B2 1525+29 






148(5) 


25/144 


3/6 


B2 2335+26 


0.25(0.05) 




6(9) 


8 


6 



Table 6. Fractions of dust morphologies. Frequency of dust morphologies as a fraction of total number of dust galaxies 
for the UGC FR-I and UGC non-FR-I sample. Columns (2)-(5) list the fractions of dust lanes, ellipses, dust structures 
which are either ellipses or lanes and irregular dust respectively. The definitions of these morphology classes are given 
in Section 0J 



Sample 


Lanes 


Ellipses 


Ambiguous 


Irregular 


UGC non-FR-I 


0.24 + 0.10 


0.28 ±0.11 


0.12 + 0.07 


0.36 + 0.12 


UGC FR-I 


0.32 + 0.13 


0.53 + 0.17 


0.10 + 0.07 


0.05 ± 0.05 



